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SUMMARY 


Model calculations of electrical currents expected to flow in Europa, which 
were summarized in a previous publication (Colburn and Reynolds, 1985) are presented 
in detail. The rationale that determined the choice of parameters is given, as well 
as the derivation of the equations defining the model, for clarification of the 
proposed current systems. 

Electrical currents should flow in the Galilean satellite, Europa, because it 
is located in Jupiter's corotating magnetosphere. The possible magnitudes of these 
currents are calculated by assuming that Europa is a differentiated body consisting 
of an outer H 2 0 layer and a silicate core. Two types of models are considered 
here: one in which the water is completely frozen and a second in which there is an 
intermediate liquid layer. For the transverse electric mode (eddy currents), the 
calculated current density in a liquid layer is approximately 10“* Am“^. For the 
transverse magnetic mode (unipolar generator), the calculated current density in the 
liquid is severely constrained by the ice layer to a range of only 10“^ to 
10"^ 1 Am - ^, for a total H 2 0 thickness of 100 km, provided that neither layer is less 
than 4 km thick. The current density is less for a completely frozen H 2 0 layer. If 
transient cracks were to appear in the ice layer, thereby exposing liquid, the 
calculated current density could rise to a range of 10” to 10 - -’ Am - ^, depending on 
layer thicknesses, which would require an exposed area of 10 - ^ to 10” of the Europa 
surface. The corresponding total current of 2.3*10^ A could in 1 yr electrolyze 
7*10^ kg of water (and more if the cells were in series), and thereby store up to 
10 1 ® J of energy, but it is not clear how electrolysis can take place in the absence 
of suitable electrodes. Electrical heating would be significant only if the ice- 
layer thickness were on the order of 1 m, such as might occur if an exposed liquid 
surface were to freeze over; the heating under this condition could hinder the 
thickening of the ice layer. 

INTRODUCTION 


Jupiter's powerful magnetic field can induce sizable electric potentials on the 
large and relatively close Galilean satellites, particularly the two innermost of 
these satellites, Io and Europa. The presence of a conducting satellite in the 
magnetosphere of its primary can give rise to two possible electric current sys- 
tems. In one of these, the transverse electric (TE) mode, eddy currents are set up 
in the satellite by the time-varying component of the magnetic field, which at 
Jupiter is caused principally by the 10° tilt of the planetary dipole field with 
respect to the rotation axis. The currents remain within the satellite body and are 
not diminished by low conductivity at the surface or outside the body. In the other 
current system, the transverse magnetic (TM) mode, a unipolar generator is set up by 
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the relative velocity of the satellite through the corotating magnetosphere (Sonett 
and Colburn, 1967). Unipolar generator currents must close within the planet's 
magnetosphere or ionosphere and consequently might be limited by the electrical 
conductivity far from the satellite. The magnetic and electrical fields appropriate 
for both the TE amd TM modes at Io and Europa are summarized in table 1 , adapted 
from Colburn (1980). 


TABLE 1.- MAGNETIC AND ELECTRIC FIELDS AT 10 AND EUROPA 


Variable 

Characteristic 

Io 

Europa 

B 0 

Magnetic field, steady component, nT 

2000 

500 

r 0 

Orbital distance, R, a 

5.91 

9.40 

v o 

Orbital velocity with respect to corotating 
magnetosphere, km sec -1 
TM mode field (E TM = v q Bq), V m“* 

56.8 

104.3 

e tm 

.114 

.0521 

V 0 

TM mode potential, per hemisphere (Vq = rE^) , kV 

207 

82 

B 1 

Magnetic field, alternating component, nT 
Angular frequency of rad sec~* 

800 

200 

ii 

0) 

1.347x10" 4 

1 .554x10 -4 

r 

Satellite radius, km 

1820 

1569 

e te 

TE mode field (E.p E = 0.5wrB.|), V m“* 

9.8lxl0~ 5 

2.44x10“-’ 


a Rj = radius of Jupiter = 7.1x10** km 


Detection of TE-mode currents requires magnetic measurements near the surface 
and measurements of or assumptions regarding the driving field. TE-mode currents 
have been measured at the Moon (Sonett et al., 1971) and used to deduce a conductiv- 
ity profile (see review by Sonett, 1982), but there has been no opportunity to 
observe these currents at Europa or Io. The TM mode is believed to have been 
observed at Io; it has been postulated as the explanation for the strong modulation 
of Jupiter's decametric radiation by Io's orbital position (Goldreich and Lynden- 
Bell, 1969), and a current was detected by the Voyager magnetometer while passing by 
Io's flux tube (Ness et al., 1979). From that experiment the magnitude of the 
current in the flux tube was estimated to be 10° A, resulting in a current of 
2x10° A at Io if northern and southern flux tubes contribute equally. Further 
analysis of the data yielded a value of 2.8x10^ A (Neubauer, 1980; Acuna et al., 
1981). For Europa, evidence for a modulation of the decametric radiation, summa- 
rized by Car and Desch (1976), is primarily negative. Among many unfruitful 
searches for the modulation effect, one ground-based study (Tiainen, 1967) showed 
increased emission when Europa was 190° from superior geocentric conjunction. 

During 1969, using RAE-1 spacecraft data, Carr and Desch found a strong Europa 
effect, also at approximately 190°, but in the following year the effect was 
indiscernible. Later Desch and Carr (1978) stated that their RAE-1 data do not show 
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a positive result concerning Europa. They cite the difficulty of sorting out Io and 
Europa effects. The approximate two-to-one ratio of the orbital periods of Europa 
and Io makes it difficult to isolate possible Europa effects in the presence of 
noise. Continuing the work of Desch, St. Cyr (1982, 1985) analyzed several ground- 
based decametric data bases, covering more than 20 yr, and found no clear evidence 
of the interaction of Europa or Ganymede with the Jovian magnetosphere. No Voyager 
measurement of flux-tube current for Europa was possible, since neither of the 
spacecraft trajectories came sufficiently close to possible flux-tube positions. 

In a detailed study, Colburn (1980) showed that internal currents from either 
mode would be very small for Io, where a conducting ionosphere envelopes a poorly 
conducting cold crust composed of such materials as sulfur, silicates, and sulfur 
compounds. Thus, nearly all of the observed current must be confined to the 
ionosphere. 

Although the induced fields are somewhat weaker at Europa (table 1), the lack 
of an observable ionosphere and the possibility of a better conducting surface 
material, 1^0 ice, overlying the liquid water, suggest that currents in the solid 
body could be larger for Europa than for Io and consequently should be investigated. 

The surface of Europa was characterized by the Voyager mission as having very 
low relief, an almost complete lack of visible craters,, many linear albedo features, 
and a surface evidently consisting almost entirely of 1^0 ice (Smith et al., 

1979). The density suggests the presence of ~6% ice by mass which, if the body is 
assumed to have differentiated, would form an outer layer of ~100 km of ice over an 
assumed iron-silicate core (Cassen et al., 1979). In this reference as well as in 
Cassen, Peale, and Reynolds (1980) the authors have further suggested the possibil- 
ity of a liquid-water layer beneath an ice crust, which could be kept from freezing 
primarily by tidally generated heat. Squyres et al. (1983) have presented a case 
for a currently active surface for Europa involving an ongoing process of cracking 
and refreezing the surface, with only a very small fraction of the satellite surface 
being fractured at any one time. The mean subsurface temperature at Europa is 
estimated to be approximately 90 K. 

In this paper we discuss electrical conductivities and calculate the TE and TM 
currents for Europa models that include a silicate interior, a surface ice layer, 
and in some cases an intermediate liquid layer. The TM-mode current densities are 
also estimated for an extreme Europa model in which surface cracks in the ice might 
allow increased current to flow. Some of the implications of such currents are 
considered. 


ELECTRICAL CONDUCTIVITY 


The electrical conductivity of rock materials has long been known to vary 
exponentially with the reciprocal of the temperature, as expressed in the relation 

o = o Q exp[-E/(nT) ] ( 1 ) 
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where o Q and E are constant over a broad temperature range, n is Boltzmann's 
constant, and T is the absolute temperature. Measurements of rock conductivity at 
a given temperature have varied by several orders of magnitude and differ apprecia- 
bly, even with different samples of the same rock type. Parkhomenko ( 1967) has 
compiled many of these determinations. More recently, the conductivity has been 
found to be influenced strongly by the oxygen fugacity and by the concentration of 
trivalent elements in the sample (Duba et al., 1974; Huebner et al., 1979). 

Although conductivity measurements can be done fairly accurately, the various 
assumptions about the material composition of a planetary body can easily produce 
uncertainties of many orders of magnitude in electrical conductivity. Colburn 

(1980) listed several rock-conductivity functions with conductivities (T □ 273) 

6 mi 

ranging from 10 to 10“ D Sm (the siemens, S, is the reciprocal ohm). The high- 
est conductivity function was ascribed to a lunar sample (Schwerer et al., 1971) and 
the lowest function was ascribed to an olivine sample (Bradley et al., 1964). 

Huebner et al. (1979) computed a conductivity function for a dunite representative 

O O 1 

of the outer lunar mantle, which gives a 273 K conductivity of 3x10 Sm . Param- 
eters for the three samples are oq = 7.9, 200, and 3.7 Sm“^ and E = 0.51, 0.92, 
and 1.09 eV, respectively. Because the computed conductivity at 273 K is essen- 
tially an extrapolation of experimental measurements, its accuracy is question- 
able. However, even over this broad range of conductivity values, it will be shown 
that currents in ice and water, when present, are always large enough that currents 
in silicates can be considered negligible in comparison. 

Measurements of the electrical conductivity of ice have varied widely, partly 
because of experimental errors. Hobbs (1974) has summarized these determinations, 
nearly all of which are taken at a temperature close to the melting point of ice 
rather than at the colder temperatures expected at the Europan surface. Measure- 
ments at nonzero frequencies (typically 1000 Hz) are not applicable to this problem, 
and dc measurements are often contaminated by potentials that develop at the elec- 
trodes. A measurement of low-temperature conductivity which appears reliable was 
made by Auvert (1973) and summarized in Auvert and Kahane (1973). Hydrogen plasma 
electrodes were used, which appears to be appropriate because laboratory experiments 
show that the method of conduction is principally proton migration (Hobbs, 1974). 
Auvert found a linear relation between the log of conductivity and the reciprocal 
temperature over the temperature range of Auvert and Kahane 's experiment, 238 K to 
158 K. For the present calculations we have used a linear fit to Auvert 's measure- 
ments and extrapolated to lower temperatures. The activation energy is 0.224 V and 
the multiplying factor is 0.0391 Sm - *'. The conductivity ranges from 3x10“^ Sm“^ at 
273 K to 1x10“ Sm -1 at 90 K. Whereas the addition of impurities to distilled 
liquid water raises the conductivity dramatically because of the addition of charge 
carriers, impurities have less of an effect on the conductivity of ice. Laboratory 
measurements of this effect are subject to error unless performed at a temperature 
cold enough to assure that no liquid phase is present. Experiments using NH^, HF, 
NHjjOH, and NH^F have been reported; use of the first three can increase ice conduc- 
tivity by one or two orders of magnitude, whereas use of NH^F decreases the conduc- 
tivity (Hobbs, 1974). It is clear that the conductivity of Europa's ice layer has 
not been satisfactorily determined. 
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If there is liquid water on Europa, it presumably contains dissolved salts from 
the silicates, as does the Earth's ocean. Salt water is much more conductive than 
distilled water. For a model of the liquid water on Europa, we use sea water of 
salinity 0.035 (standard for Earth's ocean); at 271 K (freezing point) the conduc- 
tivity is 2.75 Sm - ^ and is approximately proportional to the salinity. The conduc- 
tivity increases with temperature (Montgomery, 1963); however, the temperature 
gradient within a liquid water layer should be nearly adiabatic and hence the tem- 
perature change across the layer should be quite small. 


CURRENT DENSITIES IN THE TE MODE 


In the TE mode, the eddy currents circle the lines of the alternating magnetic 
flux (Sonett and Colburn, 1968). The geometry of the TE mode for Europa is shown in 
figure 1. The dipole comprising the principal portion of Jupiter's magnetic field 
is tilted 10° from the rotation axis; consequently a dipole component, D^., rotates 
in the equatorial plane. The time-varying field of D fc at Europa can be resolved 
into two components, in the x-direction, pointing toward Jupiter and B 2 in 
the y-direction, that of Europa 's orbital motion. The B 2 component has half the 
magnitude of B^ and lags in its time phase by one-fourth cycle, so that the 
resultant field vector traces an ellipse in Europa' s equatorial plane. An ensemble 
of eddy currents shown as I 1 encircles B-j, reaching the maximum at the time B-| 
passes through zero, and correspondingly an ensemble I 2 encircles B 2 , reaching 
the maximum at the time B 2 passes through zero. Currents 1^ and I 2 are thus 90° 
out of phase. The current density is fairly evenly distributed in any shell; the 
maximum occurs at the limb as viewed from Jupiter, whereas the minimum is less by a 
factor of only two, occurring along the sub-Jupiter meridian. For the work pre- 
sented here, the elliptical vector was approximated by a single time-varying vector 
of magnitude (B^ + B^; . 



Figure 1.- TE-mode geometry. 
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For a silicate body or one with an F^O outer layer that is completely frozen, 
the conductivities are low enough that the eddy currents near the surface are not 
modified by currents in the interior; in this case, the maximum current-density is 
given by 


i = oE ( 2 ) 

where E = 2.44x10"^ Vm" 1 and the units of i are Am" 2 . These current densities 
are plotted as a function of temperature in figure 2. The calculation is for the 
equatorial surface. For the case of salt water, if the layer were thick, the eddy 
current system would reduce the current density shown here by a modest amount. 

Calculations are shown for the three silicates previously described. For 
similar temperatures, the current density induced in ice is much larger than that 
induced in the silicates. 

If there is an intermediate liquid layer of a thickness tens of kilometers, the 
currents in the water are not significantly affected by either the silicate body 
underneath or the low- conductivity ice above. The current density for such a layer 



-J i i 

333 125 77 

T,K 


Figure 2.- Maximum TE-mode current densities as a function of surface 

material and temperature. 
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is several orders of magnitude larger than that of the solids, shown in figure 2. 

The current density in the liquid as a function of layer thickness is shown in 
figure 3, peaking at 7*10“^ Am~^ for a thin layer. The bottom of the water layer is 
at R = 1469 km. The top of the water layer is covered by ice sufficient to make 
the total radius 1569 km, but the presence of the ice does not alter the current 
density shown. The temperature of the water is assumed to be 271 K. In this range, 
the current density is almost directly proportional to the salinity of the liquid 
water. 



LAYER THICKNESS, km 


I 1 I 

1469 1519 1569 

R, km 

Figure 3.- Maximum TE-mode current density within a salt-water layer, of 
salinity 0.035, as a function of the layer thickness. 


TM-MODE CURRENT DENSITIES FOR WATER AND ICE 


The geometry for the TM mode is shown in figure 4. The steady-field component, 
B, points in the z-direction (southward) at Europa, whereas Europa's velocity, V, 
with respect to Jupiter's corotating magnetosphere is in the negative 
y(-y)-direction (i.e., opposed to the direction of Europa's orbital velocity). An 
electric field equal to V x B drives the current inside Europa in the 
-x-direction (i.e., outward from Jupiter). The direction of current flow is 
modified by variations in conductivity within Europa. If the satellite can be 
assumed to be spherically symmetric, the current will be cylindrically symmetric 
about the x-axis (symbolized by the vector I in the figure). Since Europa is 
synchronously rotating, the sub-Jupiter point is a relatively fixed position on 
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Figure 4.- TM-mode geometry. 


Europa's surface. A colatitude angle 0 can be defined as the arc between a 
point P on the surface and the sub-Jupiter point, as shown in the figure. The 
vertical current density varies as cos 0 and the horizontal current density varies 
as sin 0. The current densities discussed here are the maximum values, i.e., the 
vertical current density at 0=0 and the horizontal current density at 0 = tt/2. 

The potential which drives the TM-mode current is the product of Europa's 
diameter, the magnetic field, and the relative velocity of Europa with respect to 
the corotating Jovian-magnetic-f ield lines. The current resulting from this poten- 
tial drop must pass through Europa and then be completed by the external circuit, 
which leaves Europa in a direction nearly parallel to the field lines. Closure is 
through the distributed impedance of the Alfven wings (as discussed by Neubauer, 
1980), modified for some configurations by the Jovian ionosphere. The current in 
the interior could be bypassed by an ionosphere of Europa; however, no atmosphere or 
ionosphere in the vicinity of Europa has been reported. The current is limited by 
the total impedance. In this discussion we first consider the case in which the 
current is limited only by Europa's impedance. Then we consider the results of the 
further restriction of current by the external impedance. 

The current in the silicate part of the body is insignificant in comparison 
with possible currents in the ice or liquid layers discussed here. The satellite is 
modeled as an insulating core below 1469 km, with a liquid layer at a temperature of 
271 K topped by an ice layer of sufficient thickness to make the total radius 
1569 km. The temperature in the ice is assumed to drop linearly with distance from 
271 K to a surface temperature of 90 K, Tidal heating in the ice layer would modify 
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the shape of this profile, but not its general character. Since the ice layer 
controls the current, variations in temperature and salinity of the liquid have 
almost no effect on the current density. The resulting maximum TM-mode current 
densities are plotted in figure 5, against the water layer thickness, D w , and the 
ice layer thickness, Dj . 1^, water, is the horizontal current density in the water, 

essentially independent of depth. 1^, ice, is the maximum horizontal current den- 
sity in the ice, calculated for the interior of the ice layer, and is seen to be 
negligible in comparison to that of the water. Iy is the vertical current density 
at the ice-water interface. 

An upper limit to the total current through Europa can be estimated by compari- 
son with Io. For the Io current system, the current was apparently limited by the 
conduction of the external circuit through the flux tubes (Colburn, 1980). Neubauer 
(1980) has calculated the conduction for Europa to be less by a factor of 3 



100 50 0 

ICE LAYER THICKNESS, km 

I I I 


1469 

1519 

1569 


R, km 



Figure 5.- Maximum TM-mode current densities. 
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resulting from lower plasma densities. The driving potential for Europa is less 
than that of Io by a factor of 2. Therefore, the upper limit for the current in 
Europa can be considered to be a factor of 6 less than that observed at Io. A 
further reason for considering this to be an extreme upper limit is that the mag- 
netic field of such a current through Europa might exceed that of Jupiter and thus 
deflect the plasma and reduce the flow. We take the upper limit here to be 
2.3x10^ A, following Neubauer (1980). The corresponding maximum current density is 
3.0x10“^. Since the current densities in figure 5 are much lower than the limit 
discussed here, they are not limited by the external circuit. 


ENHANCEMENT OF CURRENT DENSITIES BY CONCENTRATION IN CRACKS 


It has been suggested that an ice layer overlaying a liquid ocean might have 
occasional vertical cracks (Cassen et al., 1979; Squyres et al., 1983). The cracks 
would be temporary because the liquid filling them would be subjected to freezing 
temperatures and to a vacuum that would make it tend to boil furiously at the sur- 
face. In the electrical context, an ice layer with an ensemble of vertical cracks 
would have an anisotropic conductivity — the vertical conductivity would be approxi- 
mated by the product of the liquid conductivity and the fractional area of the 
cracks and the horizontal conductivity would be very nearly that of the ice (assum- 
ing there are no continuous meridional paths of appreciable length on a planetary 
scale) . 


A mathematical model can be computed for this anisotropic case. The coordinate 
system chosen here is a spherical coordinate system using r, e, and <t>, with the 
applied electric field in the direction 0 = it. All of the quantities are indepen- 
dent of 4> and direction is designated by the subscripts H for horizontal and 
V for vertical; H is in the direction of +0, and V is in the direction of -r. 

The components of current density Iy and 1^ are determined by 7*1=0, 
from which 


d_ 

dr 




V 


and by v x (I/o) = 0, from which 


d_ i / d ° h\ I V g H 

dr H " ^dr J + roy “ r 


(3) 


(4) 


Since currents in the silicate portion of the body are negligible, integration 
starts at r = 1469 km with Iy = 0 and I H = 1 . After integration to 
r = 1569 km, all of the currents are scaled to satisfy v x E = 0 at the surface, 
i.e., I R = o H V 0 /r. 
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The conductivity tensor is determined from 
liquid and ice, respectively, and a factor A 
that is the open-crack area. 


o w and 


I 


the conductivities of 


surface is 


1 m 2 , 


Then for small A r 


then A c = 1 
and a 


c , which is the fraction of the surface 
As an example, if the total crack area over a 1-km 2 

• 1<T 6 . 


<< 0 , 


U 


°v = A c°w + (1 - A c )o i 


(5) 


o 


H = 


( 6 ) 


The radial variation of conductivity is approximated by 
temperature drop for the ice from 271 K at the inner surface 
surface, which determines oj. The water is assumed to have 
and temperature throughout. 


assuming a linear 
to 90 K at the outer 
the same conductivity 


By varying the fractional crack area, we calculate the horizontal current 
density in the water at the ice-water interface, I H , assuming no current limitation 
by the external circuit. The curves of I H are shown in figure 6 for three ice- 
layer thicknesses Dj = 10, 50, and 90 km. The abscissa is the proportional surface 
area of the cracks. The vertical current density at the surface Iy is also 
plotted. Values of 1^ and Iy are denoted in the figure corresponding to the 



Figure 6.- Maximum TM-mode current densities in a system with a salt-water 

layer covered by a cracked ice layer. 
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nominal upper limit for total current in Europa, 2.3x10-* A. In this case, the 
maximum Iy is 3.0x10“® Am -2 , shown on the figure as a horizontal line. The inter- 
section of this line with each curve for Iy determines the fractional crack area 
for the corresponding ice thickness. For example, the point P in the figure 
indicates that the fractional crack area for a thickness of 90 km is 6x10“^. The 
corresponding I H for this thickness is found on the I H curve (Dj^ = 90); it has 
the same abscissa as point P. Thus the three circled points on the figure show 
the I H for each thickness when the total current is at the nominal upper limit. 


The variation of maximum I H with ice thickness is better understood by means 
of an analytic approximation. Since o H of the ice is several orders of magnitude 
less than A c x o w over the range of interest, of the ice can be approximated 
as zero, whereupon the model has an analytic solution. In this solution we define 
Dj as the ice layer thickness, D w as the water layer thickness, and three radii: 
b = 1469 km, a = 1569 km, and s = 1469 km + D w . In the water (b < r < s), the 
potential is given as 


V 



(7) 


and in the ice (s < r < a), as 


V 


P + 



( 8 ) 


with constants P, Q, and M determined by the following boundary conditions: 

V = Vq (i.e. , 82 kV) at r = a and V and Iy are continuous at r = s. Then 



Ihs s2 (b 3 ' I 3 ) (* ' »)/[° w (b 3 * 

= s2lHS (?‘?)/[ a2 (?’i)] 




(9) 

( 10 ) 


or, rewriting in terms of the water and ice thicknesses D w and Dj, respectively, 


A C = I HS°' 


I v 0 - W 


v '’ [ (! : ' '■ ' (“° ' ")] 


(ID 


( 12 ) 


Here, I^g is the horizontal current density in the liquid at its upper sur- 
face, and I^Q is the vertical current density at r = a. Equation (11) has been 
written with A c on the left side for algebraic simplicity, but A c is the 
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independent variable and the equation shows implicitly the dependence of I HS on 
A q . The dependence of Iy^ on A^ follows from equation (12). The linear 
relation in figure 6 between Aq, I^g, and Iy^ is confirmed for the regime in which 
the potential across the ice (Vq - sI HS /cr w ) is essentially equal to Vq. The 
asymptotic limit for I HS suggested at the top of figure 6 is reached when the 
potential across the ice falls to zero. The near identity of the two curves for 
I H g at ice layer thicknesses of 10 km and 90 km is evident because the product 
DyDj is the same for both thicknesses, whereas s, a function of the widths, never 
varies more than 7 %. 

The maximum horizontal current density is determined by the analytic solution 
under the constraint that the total current not exceed 2.3* 10^ A; it is plotted as a 
function of liquid layer thickness in figure 7. The corresponding crack area which 
will achieve this current density is also shown. The agreement between figures 6 
and 7 is excellent for current densities as a function of liquid layer thickness, 
whereas the required crack area for any given thickness of figure 7 determined by 
the analytic solution is greater than that of figure 6 by 30# to 60#. The differ- 
ence is attributed to the neglect of ice conductivity in the analytic model. Thus, 



Figure 7.- Horizontal current density and relative crack area as a function 
of water layer thickness for the system of figure 5. 
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in figure 7, the upper Aq curve is plotted for the analytic case and the lower 
A c curve is plotted for the more realistic case in which the conductivity of ice is 
appreciable. The latter curve is constructed by translating the analytic curve 
vertically to fit the three I H points found in figure 6. Because the dependence 
of the current density 1^ on water layer thickness is nearly the same in both 
cases, it is shown as a single curve. 

-2 

It can be seen from figure 7 that 1^ will be limited to about 2x10 Am 
The relative crack area for which this can be obtained is approximately 1x10” , or 
1 m^ per 100 km^ of surface area. In this range, conduction through the ice is very 
small in comparison with that through the cracks; consequently the estimate of 
current density is not sensitive to accurate knowledge of the low-temperature con- 
ductivity of ice or the linearity of the temperature gradient through the ice. 

Since in the more general case, the abscissas in figures 6 and 7 are the product of 
the fractional crack area and liquid conductivity, any revision upward in the esti- 
mate of salinity (proportional to conductivity) would cause a lowering of the calcu- 
lated fractional-crack area by the same factor. 


ELECTROLYSIS AND HEATING FOR LIQUID-FILLED CRACKS 


The presence of Europa in Jupiter's corotating magnetosphere implies a poten- 
tial, so if the electrical conductivity is sufficient, a current should flow through 
the satellite. The magnitude of such a current would depend very strongly on the 
composition of Europa' s upper layers; the current would be very small unless a layer 
of water were present on a global scale. Among the consequences of such a current 
flow, two in particular appear to be of interest and to warrant calculation: the 

possibility of the electrolysis of H 2 0 and the amount and distribution of electri- 
cally generated heat. 


As an extreme limit, if a continuous current of 2.3x10"’ A were to electrolyze 
water in a single cell, hydrogen would be produced at the rate of 7.5x1Cr kg yr” 1 
and oxygen at the rate of 6.0x10^ kg yr” 1 . The 100-km layer of H2O contains 10 ° kg 
of hydrogen, so that significant water loss would occur only after approximately 
10^ yr. If several cells were connected in series, the rate could be correspond- 
ingly much greater; but we do not consider this possibility likely. The energy 
storage for a single cell under these conditions is lO 1 ^ J yr” 1 , while the available 
energy is 10 1 ® J yr -1 . 


The possibility of electrolysis depends on the electrical-conduction mecha- 
nism. Laboratory experiments show that conduction in ice occurs principally through 
the migration of protons. In one experiment, when the electrodes were made of a 
porous black metal saturated with hydrogen, protons were provided at the anode by 
ionization of hydrogen, whereas at the cathode, protons were neutralized to produce 
hydrogen. The current decreased over the course of time as th4 anode became more 
depleted, and the experimenter postulated that eventually the protons would be 
supplied by natural ionization of H 2 0 (Hobbs, 1974). Thus, oxygen production by 
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electrolysis would not be favored in a hydrogen environment. Auvert's experiment 
used hydrogen plasma for both electrodes, and consequently protons were readily 
available to enter the ice at the anode and leave at the cathode. The external 
circuit provided the electron transfer to maintain charge neutrality. 

At Europa the interfaces considered likely are (1) between the surface of 
Europa and the tenuous atmosphere and (2) between ice and seawater. The atmosphere 
is ionized gas, which is electrically conducting, consisting largely of protons and 
electrons, so that hydrogen plasma electrodes similar to Auvert's could be main- 
tained. The protons supplied by the plasma could carry the charge through the ice 
and cause proton depletion at the anode and concentration at the cathode, while the 
external circuit could provide the electrons needed for charge neutrality. However, 
appreciable concentrations of plasma would not be expected to occur since the plasma 
corotates with Jupiter's field rather than being tied to Europa, and general turbu- 
lence and thermal effects would maintain a relatively uniform plasma concentration 
over time, providing a reservoir of protons. Thus, we cannot expect either hydrogen 
or oxygen to accumulate at the icy surface. If there were liquid water at the 
surface, chemical interactions between the solvated ions and plasma could possibly 
allow hydrogen and oxygen generation, but this interface would last only a short 
time before freezing over at the prevailing low surface temperature. 


At the water-ice boundaries, it is questionable whether gaseous hydrogen or 
oxygen would be generated. The protons leaving the ice could dissolve in the water 
as H + ions. At the other boundary, where protons need to enter the ice, H + ions 
(which are always present in water) could leave the water. Thus, in the absence of 
metallic electrodes, electrolysis would probably not occur. 

Drobyshevski (1980, 1982) has proposed a method called "volume electrolysis" 
whereby unipolar generator currents in icy planets such as Europa could produce 
large quantities of hydrogen and oxygen in the interior, which could subsequently be 
ignited upon meteoric impact. He suggested that over time, the buildup of the gases 
could be large enough that the explosions could be of near-planetary proportions and 
remove sufficient material to change the satellite mass significantly. If the 
current were limited to a value commensurate with that observed at Io (Drobyshevski 
used 10? A), the amount of gas buildup discussed in his model could occur only if 
the current were to pass through numerous cells in series to multiply the gas pro- 
duction. Drobyshevski proposes that inhomogeneities in the ice would provide these 
cells; although it is not obvious how such a series of cells could be produced, the 
current could pass around such cells rather than being forced through them. In any 
event, if the electric field were similar to that now present at Europa, ohmic 
losses in the ice should reduce the current by many orders of magnitude from that 
which he proposed. 


The total available electrical power averaged over the surface is only 
10 -3 Wm - ^, which is insignificant when compared with the 50 Wm - ^ of solar energy 
arriving at the surface. In a discussion of possible electrical currents in Io, 
Gold (1979) suggested a positive-feedback mechanism in which any tendency for cur- 
rents to concentrate in a local area would be enhanced by localized electrical 
heating which would reduce the resistivity and, in turn, draw more current. Under 
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proper conditions, the same process might occur at Europa. However, this possibil- 
ity must be weighed against the extra heat losses that such a concentration would 
cause. We consider first whether such a concentration could maintain itself if 
created by fissures in the ice. If a crack occurred that exposed liquid at Europa's 
surface, some of the energy removed by boiling might be provided by electrical 
heating, although heat losses would be so great that an open, liquid surface would 
tend to freeze rapidly. When water boils in a vacuum at 273 K, the flux of escaping 
water is about 0,5 kg m - ^ sec - \ as determined from a one-dimensional momentum- 
transport calculation. 


The details of this calculation are as follows. In the absence of another gas, 
water will boil at 273 K if the pressure of the vapor above the liquid surface is 
610 kg m - ^ sec - ^, as given in vapor pressure tables. The average speed of water 
molecules at this temperature is 615 m sec"^ (since the velocity (v) is defined by 
v^ o 3 k T/m, where k is Boltzmann's constant and m the mass of a molecule). 

The momentum blocked by a wall of the vapor container is half the pressure, or 
305 kg m~^ sec”^. Thus, if the wall were removed, a rough approximation is that a 
flux of 0.5 kg m - ^ sec - would be maintained by the boiling process. The back 
pressure of the vapor would reduce this flux, but the number here represents a first 
order approximation. 


The heat of vaporization for this flux is 1.3x10** W m - ^, 


and the power lost by 
Again, for a current comparable to that at 


atellite is 4x10^® W, which is enough to 
m^ of exposed water surface under the most 


radiation is negligible in comparison. 

Io, the available power for the whole 
maintain by electrical heating only 10 
optimum conditions. For such a concentration of current, the plasma surrounding 
Europa would have to be able to conduct the current to the cracks with low loss, 
which would require unusually high plasma conductivity. Hence, it appears that an 
exposed liquid surface could not be maintained by electrical heating. 


A situation in which electrical heating might impede the closure of a crack 
involves the conditions just after an initial ice layer has formed, with a column 
mass sufficient to contain the vapor pressure of the water, i.e., a thickness of 
approximately 0.5 m. If the radiative heat loss to space and the heat loss by 
horizontal conduction could be provided electrically, further freezing would not 
occur. This condition would require that the ice remain at a relatively high tem- 
perature (near 273 K), for if it were to cool significantly, the electrical conduc- 
tivity would become very low and would shut off the available power. The heat lost 
by radiation from a 273-K ice surface is ~300 Wm - ^. With the ice-column resistance 
of 1.7x10-* ohms m - ^, this heating rate could be produced by a current of 0.04 Am - ^, 
thus maintaining a thin crack-surface of 5 km^ at constant thickness. The voltage 
drop across such a plug would be only 10% of the total available for this case, so 
that salt-water conductivity would not limit current flow. The heat loss to the 
surrounding ice would be less than the surface loss for horizontal thermal gra- 
dients, in which the entire temperature drop of about 1 80 K occurs over a distance 
of a few meters or more. Such cracks would also have to be maintained simulta- 
neously in both hemispheres if a current were to flow. It is assumed here that the 
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cracks had been created by tidal action or by the pinching effect described by Gold 
(1979) (if such a mechanism can operate at the temperatures of Europa's surface). 


CONCLUSIONS 


Electrical currents through Europa would be expected to be orders of magnitude 
larger than those through Io because of the lack of a significant ionosphere and the 
relatively high conductivity of ice in comparison with silicates. However, the 
currents are still shown to be quite small for a solid satellite. If a liquid water 
layer were to exist on Europa of a salinity comparable to that of the Earth's 
oceans, a current density of 10“^ Am -2 should flow in the TE mode, since the flow 
depends only on the strength of the alternating magnetic field and the conductivity 
of the water, both of which are known. However, electrical heating would be negli- 
gible at this current density level and electrolysis would not occur because the 
complete circuit is in the liquid which contains no electrodes. On the other hand, 
the current density in the TM mode depends on a larger number of factors of which 
only the E = V x B potential is known. The magnitude of the current density 
depends on the conductivity in the total circuit comprising the liquid, the ice, the 
flux tubes in the magnetosphere, and a connection somewhere in the neighborhood of 
Jupiter's ionosphere; but this conductivity could be much less than in the corre- 
sponding circuit at Io. It is thus quite possible that no significant TM-mode 
currents can occur. It is shown for an ice layer at least 4 km thick, the current 
density in the liquid will be restricted to 10" ^ Am -2 or less. It is also shown 
that if cracks occurred in the ice layer that exposed 10~9 to 10”® of the Europa 
surface to the liquid layer, average current densities as high as 2x10“^ Am -2 could 
be reached. If electrolysis were to occur under these conditions, with a total 
current of 2.3*10^ A, then 7x10^ kg yr - ^ of water could be electrolyzed with an 
annual energy storage of 10 ^ j. However, it is not certain that the interfaces 
between the liquid, ice, and external plasma provide suitable electrodes for elec- 
trolysis. The electrical heating from TM-current levels discussed here is inconse- 
quential in comparison with other energy sources on a planetary scale, and is shown 
to be significant only when concentrated at an ice layer approximately 1 m thick, 
such as might occur after the freezing over of a crack in an ice layer. The Galileo 
spacecraft mission to Jupiter will make several close fly-bys of Europa and should 
permit the measurement of any significant current through the Europa flux tubes. 
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